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Abstract 

Increases in drought and temperature stress in forest and woodland ecosystems 
are thought to be responsible for the rise in episodic mortality events observed 
globally. However, key climatic drivers common to mortality events and the 
impacts of future extreme droughts on tree survival have not been evaluated. 
Here, we characterize climatic drivers associated with documented tree die-off 
events across Australia using standardized climatic indices to represent the key 
dimensions of drought stress for a range of vegetation types. We identify a 
common probabilistic threshold associated wdth an increased risk of die-off 
across all the sites that we examined. We show that observed die-off events 
occur when water deficits and maximum temperatures are high and exist out- 
side 98% of the observed range in drought intensity; this threshold was evident 
at all sites regardless of vegetation type and climate. The observed die-off events 
also coincided with at least one heat wave (three consecutive days above the 
90th percentile for maximum temperature), emphasizing a pivotal role of heat 
stress in amplifying tree die-off and mortality processes. The joint drought 
intensity and maximum temperature distributions were modeled for each site 
to describe the co-occurrence of both hot and dry conditions and evaluate 
future shifts in climatic thresholds associated with the die-off events. Under a 
relatively dry and moderate warming scenario, the frequency of droughts capa- 
ble of inducing significant tree die-off across Australia could increase from 1 in 
24 years to 1 in 15 years by 2050, accompanied by a doubling in the occurrence 
of associated heat waves. By defining commonalities in drought conditions 
capable of inducing tree die-off, we show a strong interactive effect of water 
and high temperature stress and provide a consistent approach for assessing 
changes in the exposure of ecosystems to extreme drought events. 



Introduction 

Drought is a pervasive feature of forest ecosystems that 
constrains primary productivity (Zhao and Running 
2010) and, during extreme conditions, can induce large- 
scale dieback (loss of above-ground tissues) or mortality 
episodes (Breshears et al. 2005). A global surge in the 
study of drought-related impacts on forests has helped to 
document an increasing number of mortality events in 
the literature and implicate temperature increases as 
amplifying moisture deficit, heat stress, and the impacts 
of biotic agents on tree species (AUen et al. 2010; Toomey 
et al. 2011; Williams et al. 2013). In Australia, forest and 



woodland ecosystems are strongly influenced by large cli- 
matic variability and recurring drought events. Recently, 
these drought patterns have been affected by systematic 
shifts in precipitation and rising temperature. Indeed, 
recent drought-related forest die-off events observed in 
southwestern Australia (Matusick et al. 2013) were 
accompanied by increases in mean temperatures and the 
incidence of short periods (>2 days) of temperature 
extremes or heat waves likely to induce acute heat stress 
preceded by a long-term (~40 years) decline in mean 
annual rainfall (Fig. SI). Drought-induced forest die-off 
could either represent episodic disturbances within an 
existing climate regime or be indicative of a climate shift. 
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While episodic disturbances are likely to facilitate stability 
of forest structure and function over the long-term 
(Lloret et al. 2012), climatic shifts wiU induce significant 
changes in existing forest dynamics and the potential for 
state changes (Rietkerk et al. 2004) The challenge in dis- 
entangling which of these two processes dominates lies in 
defining the frequency, intensity, duration, and trajectory 
of climate extremes capable of inducing die-off so that 
future scenarios may be evaluated based on the likelihood 
of drought events. 

Drought is defined as a prolonged and exceptional defi- 
cit between water supply and demand, in our case at the 
tree or stand level. It can be characterized by both its 
duration and intensity and is often associated with periods 
of above average temperatures and heat waves (Dery and 
Wood 2005). Species and communities may be differen- 
tially impacted by drought (Fensham and Fairfax 2007; 
Koepke et al. 2010; Mitchell et al. 2008), being dependent 
on the primary meteorological drivers, and secondary fac- 
tors such as the presence of biotic agents and conditioning 
factors that determine site and species sensitivity to 
drought (Mitchell et al. 2013a). The mechanisms by which 
drought intensity and duration mediate species sensitivity 
through the physiological processes controlling plant 
carbon and water balance remain unresolved (McDowell 
et al. 2008; Sala 2009). Recent conceptual models have 
implicated two interrelated physiological mechanisms: 
hydraulic failure (desiccation of water conducting tissues 
within the plant) and carbon starvation (depletion of 
available carbohydrates and failure to maintain defenses 
against biotic agents) (McDowell et al. 2008, 2011) as 
being the primary mechanisms through which trees may 
succumb to drought stress. Within this framework, the 
relative importance of these two underlying physiological 
processes is intimately linked to the key attributes of the 
drought regime; intensity, duration, and frequency and 
how biotic and abiotic factors modulate plant water avail- 
ability through time (Anderegg et al. in press). 

Rising temperatures and the associated increase in the 
frequency of heat waves further increases the risk of mor- 
tality by altering plant water and carbon balances or by 
increasing the presence and activity of biotic agents, 
regardless of changes in precipitation regimes. Increasing 
evaporative demand will intensify drought conditions in 
the root zone and/or plant canopy and accelerate dehy- 
dration resulting from irreversible embolism of the xylem 
(Brodribb and Cochard 2009; Sperry et al. 2002). 
Increases in plant growth, tree size and possibly leaf area 
index associated with elevated temperature and CO2 con- 
centrations under favorable conditions could intensify 
intertree competition (Warren et al. 2011) and increase 
forest vulnerability to mortality during sudden or pro- 
tracted periods of drought (Duan et al. 2013; Levanic 



et al. 2011). Rising temperatures can also increase the 
activity of biotic agents by weakening plant defenses 
(Boyer 1995) or provide suitable conditions for new pest 
species (Kurz et al. 2008), thereby increasing the severity 
of plant stress during drought. The impacts of elevated 
temperatures on plant functioning represent a critical 
dimension of drought risk, and the potential for interac- 
tive effects operating in concert with the consequences of 
water deficit therefore need to be considered when char- 
acterizing drought. 

The role of climate extremes may be critical in shaping 
future species and ecosystem dynamics (Zimmermann 
et al. 2009), but because of the sporadic and unpredict- 
able nature of these extreme events, they remain difficult 
to document and monitor through time and space. The 
manner in which drought is defined and described in 
relation to ecosystem processes tends to be context spe- 
cific (Smith 2011) and restricts our ability to compare the 
climatic conditions associated with reported episodic die- 
off events from different regions and climates. Recent 
evidence from global meta-analyses of hydraulic safety 
among woody species from a range of climates suggest 
that many species' water transport systems have evolved 
to operate very close to the limits imposed by their envi- 
ronment (Choat et al. 2012). This physiological evidence 
and the occurrence of episodic forest die-off events across 
a wide range of forest types (Allen and Breshears 1998; 
Keith et al. 2012) would suggest that climate thresholds 
on plant processes such as mortality may be reached at 
similar extremes relative to the long-term climatic distri- 
bution of relevant climatic parameters. 

The documented occurrences of tree die-off across dif- 
ferent ecosystems provide a priori evidence that critical 
thresholds have been surpassed with respect to plant 
functioning and health. However, the lack of consistent 
ecological information available from historic events limits 
the applicability of predictive modeling techniques in 
many cases. An alternative approach is to use the primary 
drivers, that is, the climatology to build a consistent pro- 
file of the drought responsible for inducing known occur- 
rences of tree die-off. This is akin to environmental 
envelope techniques used to describe presence-only 
records of species (Pearce and Boyce 2006). We sought to 
develop a suitable approach that represents key climate 
drivers of these events within a probabilistic framework 
such so that drought attributes are described in terms of 
deviations from long-term mean conditions irrespective 
of underlying differences in climatology. The approach 
also defines the mean intensity, duration, and frequency 
of discrete events given their relevance to plant function- 
ing and the underlying mechanisms responsible for tree 
mortality. This approach allowed us to address the 
following questions: (1) What are the key climatic 
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thresholds of documented events based on commonalities 
in relevant indices among different sites? (2) What is the 
extent to which high temperature co-occurs with high 
levels of water deficit? (3) What are the projected changes 
in the frequency, intensity, and duration of extreme 
drought based on a relatively dry and moderate warming 
scenario for Australia? 

Materials and Methods 

Australia is a predominately water-limited environment, 
and ecosystem productivity is tightly coupled to precipita- 
tion (Adams 1995). Here, we used Australia as a case 
study to investigate the predisposing climatic conditions 
associated with documented die-off events identified from 
a survey of published data. The majority of global circula- 
tion models (GCM) predict declines in precipitation 
across much of Australia coupled with increases in mean 
annual temperatures and exceptionally hot periods (Hen- 
nessy et al. 2008; Lucas et al. 2007). Our probabilistic 
approach uses the standardized precipitation evapotrans- 
piration index (SPEI) and a similar maximum tempera- 
ture index (MTI; see below) to define those drought and 
temperature stress thresholds (based on univariate and 
joint distributions) and drought attributes (intensity, 
duration, and frequency) associated with historic die-off 
events across a range of forest and woodland communi- 
ties. We then evaluate the likelihood of future conditions 
surpassing these thresholds using a moderate warming 
scenario with GCM-derived climate data (2011-2050) and 
finally conduct a sensitivity analysis based on observed 
climate (1961-2010) scaled to match future patterns in 
temperature and precipitation at 2050. 

Drought die-off events 

We conducted a literature search using various online 
databases to identify incidences of tree mortality and die- 
back that were primarily attributed to drought or where 
there was evidence that drought was the trigger for tree 
death or canopy collapse. Die-off events were identified 
from on-ground measurements of stand health that 
showed evidence of mortality or canopy collapse across a 
range of size classes, foDowed by significant increases in 
mortality during or directly after the event. For all events, 
drought damage or mortality was estimated at >5% of 
the individuals in relatively mature-aged stands, using a 
range of sampling strategies. Two sites (Canberra and 
Hobart) had multiple drought events, and their thresh- 
olds were assessed for each event, that is, 15 sites and 17 
events in total. Data on the severity and extent of the 
reported tree die-offs ranged from long-term plot inven- 
tory measurements that had been repeatedly sampled 



through time (Cunningham and Walker 1973) to oppor- 
tunistic observations of die-off (Took et al. 1966). The 15 
documented sites included in the analyses ranged fi-om 
wet sclerophyll forest, to open woodland and shrub lands 
that were generally dominated by Eucalyptus spp. or Aca- 
cia spp. (mean annual precipitation 240-1161 mm) and 
covered many different biomes across Australia (Table 
S2). Biotic agents such as wood-boring insects were found 
to be present in about half of the die-off events reported 
(Table S2). 

Climate analysis 

Daily meteorological data were extracted from the meteo- 
rological station closest to each of the 15 sites (<30 km 
away) using patched point datasets - observed station 
data where missing or suspect values are "patched" with 
interpolated data (http://www.longpaddock.qld.gov.au/ 
silo/) (Jeffrey et al. 2001). Among all 15 stations, an aver- 
age of 20% and 50% of the rainfall and temperature data, 
respectively, were "patched" following the methods 
outlined by Jeffrey et al. (2001). Monthly climatic water 
deficits (1891-2012) were calculated using the SPEI (Vi- 
cente-Serrano et al. 2009). The SPEI is a simple climatic 
water balance based on the difference between precipita- 
tion and potential evapotranspiration (calculated from 
Thornthwaite's equation), that is, de-seasonalized and 
subsequently standardized, so that drought conditions at 
different locations, over a variety of predefined time 
scales, can be represented as probabilities and presented 
in a comparable fashion. The formulation of SPEI helps 
to remove the seasonal autocorrelation in the water bal- 
ance time series by fitting a three parameter log-logistic 
distribution separately to each calendar month. The stan- 
dardization technique allows for direct comparison of 
SPEI values among sites and transforms the fitted data to 
standard normal, that is, mean is approximately zero and 
the standard deviation is one (Vicente-Serrano et al. 
2009). Rather than adopting the full physical formulation 
of evapotranspiration provided by models such as Pen- 
man-Monteith, the Thornthwaite formulation of potential 
evapotranspiration was used to estimate monthly climatic 
water deficit due to inherent limitations in existing and 
future datasets (particularly the availability of daily relative 
humidity and net radiation data at all sites). We tested 
whether the Penman-Monteith equation (FAO 56, Allen 
et al. 1998) produced significantly different patterns in 
SPEI among the 15 sites over the observed 121 year time 
series. The SPEI values from the two methods were highly 
correlated (R^ = 0.95-0.99), and their slopes were not sig- 
nificantly different from one. Furthermore, both simula- 
tions exhibited similar minimum SPEI values during the 
documented die-off events. We used a SPEI calculated for 
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a 6-month analysis window for the identification and 
computation of drought events and their attributes using 
the methods described in Vincente-Serrano et al, (2009). 
The 24 and 48 months SPEI values were also calculated to 
assess whether there was a longer-term signal in drought 
conditions prevailing at the start of the documented die- 
off events. By construction, the SPEI has a standard nor- 
mal distribution; hence, negative values represent periods 
of net negative water balance, and the probability of 
exceeding values of —1 and —2 is approximately 0.16 and 
0.02, respectively. 

We developed a maximum temperature index (MTI) 
based on the monthly 90th percentile of daily maximum 
temperature. The MTI provides a monthly estimate of high 
temperature stress and is similarly standardized to allow 
for comparison across different climates. Specifically, we 
recorded the daily maximum air temperature from the 
observed climatic record (1891-2012), extracting from this 
the 90th percentile for each month. We de-seasonalized 
this by fitting a linear model with sinusoidal covariates on 
a 12-monthly period. After checking for trend, we fitted a 
three parameter log-logistic distribution to the residuals 
using the method of L moments, double-checked parame- 
ter estimates with maximum-likelihood estimates, and 
subsequently transformed the probability distribution of 
the log-logistic quantUes to a standard normal distribution. 
The process of deseasonalizing the data, fitting log-logistic 
distributions to the residuals, and then transforming these 
to standard normal was identical to the procedure used in 
the construction of the SPEI (Vicente-Serrano et al. 2009). 

The joint distributions of SPEI and MTI were fitted 
using three copulas: Gaussian, Student-f and Frank. Cop- 
ulas were used because they are capable of capturing 
linear and nonlinear forms of dependence between their 
marginal inputs and are therefore a more powerful tool 
for dependence analysis than say correlation coefficients. 
A Gaussian copula fitted to marginal Gaussian distribu- 
tions (as is the case here) returns the well-known multi- 
variate Gaussian distribution, characterized by the linear 
correlation matrix Q. We used this to represent the "null 
hypothesis" of simple linear dependence between SPEI 
and MTI. The Student-f copula is parameterized by the 
linear correlation matrix Cl and a degrees-of-freedom 
parameter v. This copula, like the Gaussian, is symmetri- 
cal, and the Gaussian copula is a special case of the 
Student-f as v approaches infinity. Importantly, however, 
at low degrees of freedom the Student-f copula allows for 
tail dependence - that is a greater degree of dependence 
between the marginal covariates as their values become 
much higher, or much lower, than their means. The 
potential for tail dependence in this context is important 
because this would indicate a greater probability of 
vegetation experiencing both high water deficit and high 



temperature stress than what would occur under a normal 
(Gaussian) distribution. The Frank copula is also parame- 
terized by a single dependence parameter and serves as an 
alternative to the Gaussian, particularly as it lacks tail 
dependence and when fitted to standard normal margi- 
nals displays greater dependence around the marginal 
means. In other words, this implies that the probability of 
conditions of high water deficit and temperatures are less 
than what would be expected under a normal distribu- 
tion. All copulas were fitted using the "copula" library in 
the R program (v2.11.1, R Foundation for Statistical 
Computing, Vienna, Austria) with maximum-likelihood 
parameter estimation. The final copula for each site was 
chosen using a simple likelihood ratio test (see Table S3). 

Drought events were identified from the time series of 
SPEI using a 6-month timescale. We used this timescale 
because meteorological drought can be as much as 
6 months out of phase with soil water deficit (Mpelasoka 
and Chiew 2009) and droughts shorter than 6 months are 
unlikely to place severe stress on the predominately scle- 
rophyllous vegetation of Australia. A drought event was 
defined as occurring whenever the SPEI values were <— 1 
following (McKee et al. 1993). The drought commenced 
when the SPEI value became less than zero and ended 
when the values had recovered back to zero or higher. 
Drought duration was calculated by summing the months 
during this period, and drought intensity was taken as the 
mean SPEI over the drought period (see Fig. SI for an 
example). The relationship between antecedent conditions 
at the beginning of each die-off event (SPEI at 24 and 
48 month intervals) and its duration and intensity was 
used to test whether the previous long-term water balance 
had any effect on the drought attributes during the event 
using linear regression (at P < 0.05). Heat waves were 
identified as three consecutive days between October and 
March over the 90th percentile threshold (based on a 15- 
day analysis window) and were calculated using the 
"CTX90pct" method outlined in Perkins and Alexander 
(2012). In their appraisal of a large range of heat-related 
indices, Perkins and Alexander (2012) showed that the 
CTX90pct method to calculate heat waves was suitable 
for capturing trends in heat waves in Australia and pro- 
duces an adequate population of measurable events that 
are optimal in describing "extreme" versus "measurable" 
periods of heat stress. 

Climate projections were generated for the 15 sites 
using the CSIRO Mk 3.5 GCM and A2 emissions scenario 
at a medium sensitivity (WCRP CMIP 3, http://esg.llnl. 
gov:8443/data/). Intermodel comparisons show that this 
model produces moderate warming and some of the larg- 
est declines in annual rainfall across Australia and 
performs very well for the key features of the Australian 
climate (Watterson 2012) (see Fig S4). The requirements 
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of daily GCM data over the analysis period restricted the 
choice of suitable models for the analysis. Daily data of 
precipitation and maximum and minimum surface tem- 
perature were downscaled to the fifteen sites using two 
different approaches. Firstly, historic GCM precipitation 
data (1891-2010) from the nearest grid point were 
matched to observed precipitation data (1891-2010) using 
a daUy translation method (Mpelasoka and Chiew 2009). 
Secondly, daily temperatures from the GCM grid cell were 
scaled to the observed station data using a constant factor 
(ratio of daily mean GCM historic to observed tempera- 
tures) and then applied to the GCM future data. SPEI 
values, MTI values, and drought event characteristics 
could then be computed for the period between 1891- 
2010 (historic) and 1891-2050 (future), based on a con- 
tinuous monthly time series of precipitation and potential 
evaporation. A sensitivity analysis was conducted to test 
for the effect of projected changes in temperature and 
precipitation independently. Projected changes in temper- 
ature and precipitation were applied to the observed sta- 
tion data (1961-2010) using monthly scaling factors for 
maximum and minimum temperature (°C change per °C 
of global warming) and precipitation (percentage change 
per °C of global warming). All analyses were performed 
using the R program and relevant packages (v2.11.1, R 
Foundation for Statistical Computing). 

Results 

Using the observed climate record (1891-2012), our 
analysis showed that 17 die-off events across a range of 
vegetation types occurred between 1941 and 2012 when 
conditions were both exceptionally dry and hot (Fig. lA- 
G, Table 1). The impact (% of stand affected ranged from 
-7-90%), total extent, and spatial patterning of these 
events were not consistent across the landscape with die- 
off occurring along guUies and riparian areas as well as 
shallow, exposed slopes and hill tops (Table S2). During 
die-off events, drought duration ranged from 8 to 
50 months and mean intensity (mean SPEI) ranged from 
— 1.00 to —1.54 (more negative value represents higher 
drought intensity), highlighting the range, duration and 
intensity of conditions capable of inducing die-off 
(Fig. 2A, Table 1). Among all observed die-off events, the 
monthly sequence of SPEI values reached probabilities of 
<0.02 for 1 month or more, regardless of whether the 
event was of high intensity/low duration or low intensity/ 
high duration (Figs. 1 and 2 A, Table 1), suggesting that 
extreme conditions of this magnitude or greater were 
required to induce drought die-off 

For 14 of the 17 die-off events, the MTI values reached 
probabilities > 0.98 of the observed range during the 
drought event (Table 1), while all events surpassed MTI 



probabilities > 0.97. At least one heat wave (see Methods 
for definition) occurred during each of the documented 
tree die-off events (Table 1). Among all drought events in 
the sites' climate records, 99% of "extreme" droughts (see 
definition below) contained at least one heat wave high- 
lighting the relevance of short stochastic periods of tem- 
perature stress to vegetation experiencing water deficit. By 
analyzing the conditions common to all die-off events, we 
classified "extreme" droughts as those that included 
months that were <0.02 percentile for SPEI, >0.98 percen- 
tile for MTI and included at least on heat wave event. 
This definition allows us to distinguish between more 
common drought events (where SPEI < 0.33 percentile) 
that occur, on average, 1 in 3 year, from extreme 
droughts that occur less frequently at 1 in 24 year for the 
observed climates among the 15 documented sites 
(Table 2). 

The influence of the SPEI sequence prior to docu- 
mented die-off events was examined to assess whether 
die-off events were consistently associated with a particu- 
lar set of antecedent conditions (i.e., exceedingly dry or 
wet conditions). Documented die-off events occurred 
after both wet and dry periods; in all observed die-off 
events the 24-month SPEI (at the start of the drought 
event) ranged from —2.26 to 1.57 (Table 1). Furthermore, 
the duration and intensity of extreme drought events 
within the climate record were never correlated with 
either the 24- or 48-month SPEI based on our regression 
analyses (data not shown). 

The mean monthly SPEI versus MTI values during die- 
off events were plotted on the fitted joint probability den- 
sity function of observed SPEI and MTI (Fig. lA-G). For 
about one-third of the sites, a likelihood ratio test sug- 
gests that the best fit of the joint distributions for the 
observed climates were Student- f or Frank copula distri- 
butions. For these sites, the fitted copula distributions 
indicated strong ( Student- f) or weak (Frank) tail depen- 
dence between the SPEI and MTI values (Fig. 1B,C and 
E, Table 3). The joint probability density function using 
precipitation and temperature predictions from a moder- 
ate warming and dry scenario were also fitted to show 
changes in drought conditions (Fig. 1, Table 3). Future 
climates at aU sites showed a trend toward greater proba- 
bility of hotter and drier drought conditions, and the 
probability of sites experiencing the threshold correspond- 
ing to extreme conditions (defined as SPEI < 0.02 percen- 
tile and MTI > 0.98 percentile) increased by 56-270% 
across the different sites (Table 3). 

Under the relatively moderate warming scenario with 
associated declines in precipitation used here, drought 
events became longer, more intense and/or more frequent 
(Fig. 2A-D, Table 2). In general, the frequency of pre- 
dicted drought events relative to the historic data 
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Figure 1. Probability of observed and future changes in drought intensity and maximum temperature. (A-G) Subset of tree die-off sites from a 
range of vegetation types across Australia showing fitted joint probability density of SPEI and Mil for observed (1891-2010; shaded background 
and dashed lines) and future projections (2011-2050 based on CSIRO Mk 3.5, A2 SRES scenario; solid contour lines). The three contour lines 
denote probability densities at 0.15,0.05, and 0.02. The mean SPEI and MTI value (green square) and corresponding maxima and minima (black 
crosses) during the documented drought-induced die-off events at each site is plotted (see Fig. SI for an example of how these were derived). 
Dashed vertical and horizontal lines indicate probabilities of SPEI < 0.02 percentiles and MTI > 0.98 percentiles (for their unimodal distributions) 
under observed climate conditions. The copula distribution fitted to the observed and future projection is given in the lower right for each panel. 
The top left hand inset (H) gives an overview for interpreting the joint probability density of monthly drought intensity and maximum temperature 
at any individual site. The trajectory of future climate is shown by arrows indicating that climate may change in terms of both temperature and/or 
drought (combination of precipitation and temperature). (I) Map of Australia showing the distribution of projected percentage changes in 
precipitation for 2050 (using CSIRO Mk 3.5, A2 SRES scenario) based on an observed baseline (1975-2004) and locations of die-off events 
denoted by lower case letters corresponding to panels A-G. 



increased by an average of 28% (—12 to 79%) across all 
sites with a 20% and 7% increase in duration and inten- 
sity respectively (Table 2). More importantly, the fre- 
quency of events considered extreme and capable of 
inducing tree die-off increased at all sites by 65% (—36 to 



234%), changing from occurring at 1 in 24 years to 1 in 
15 years (Table 2). The incidence of heat waves also dou- 
bled and droughts were more likely to co-occur with heat 
waves (Table 2). A small number of predicted droughts 
fell outside the margins of the observed record and are 
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Table 1. Details of the drought attributes associated with documented die-off events in Australia. Drought intensity represents the mean of all 
monthly SPEI values during the event. The monthly minimum SPEI and maximum MTI values represent extreme water deficit and high temperature 
conditions during the event, and their corresponding percentiles are given in brackets. Superscript letters after the site name refers to panels 
shown in Fig. 1 . 







Duration 


Intensity 


Min SPEI 


Max MTI 


No. of 


Predrought 


Site 


Start year 


(months) 


(mean SPEI) 


(min percentile) 


(max percentile) 


heat waves 


SPEI (24 month) 


Alpha'' 


2001 


46 


-1.16 


-2.57 (0.007) 


2.66 (0.999) 


17 


-0.55 


Armidale"^ 


1982 


16 


-1.47 


-2.46 (0.003) 


2.29 (0.997) 


16 


-2.26 


Bollon 


1979 


29 


-1.02 


-1.86 (0.020) 


2.04 (0.993) 


14 


0.76 


Canberra'' 


1965 


8 


-1.38 


-2.24 (0.005) 


2.13 (0.987) 


4 


0.59 


Canberra'' 


1981 


20 


-1.26 


-2.05 (0.011) 


2.37 (0.998) 


16 


-0.62 


Charters Towers" 


1991 


50 


-1.16 


-2.26 (0.003) 


2.23 (0.995) 


25 


1.57 


Cobar 


1965 


18 


-1.18 


-1.82 (0.019) 


2.35 (0.997) 


4 


0.85 


Cooma 


1965 


8 


-1.09 


-1.83 (0.016) 


2.13 (0.988) 


4 


0.95 


Hobart* 


1977 


47 


-1.09 


-2.27 (0.006) 


1.97 (0.988) 
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15 
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1.74 (0.971) 


8 


0.99 
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2010 


22 


-1.50 


-2.70 (0.001) 


2.10 (0.995) 


6 


1.27 


Mt Macedon'^ 


1967 


14 


-1.24 


-2.09 (0.010) 


2.11 (0.989) 


7 


0.05 
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1967 


8 


-1.51 


-3.21 (0.001) 


2.37 (0.991) 


1 


0.44 


Yeelirrie 


1976 


21 


-1.00 


-1.98 (0.020) 


1.74 (0.969) 


5 


0.91 
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2002 


13 


-1.40 


-2.14 (0.004) 


1.80 (0.976) 


7 


0.99 


Wilcannia 


1941 


11 


-1.09 


-2.04 (0.017) 


1.96 (0.983) 


4 


0.20 




Mean 


22 


1.25 


-2.22 (0.009) 


2.10 (0.989) 


9 


0.44 
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8 


-1.54 
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1.39 (0.969) 


1 


2.26 
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50 


-1.00 


1.82 (0.020) 


2.66 (0.999) 


25 
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perhaps indicative of "mega-drought" conditions, charac- 
terized by higher intensities and longer durations than 
has been observed in the historic record (Fig. 2A). 

Given that predictions of precipitation are particularly 
variable across GCMs and scenarios, we performed a sen- 
sitivity analysis to help distinguish the role of temperature 
and precipitation on drought attributes. Simulated 
changes to observed records (1961-2010) of temperature 
(-1.44 °C change from 2010 climate) and precipitation 
(11-27% annual decline by 2050) were found to increase 
drought event duration by approximately 22% with little 
change in mean drought intensity (Fig. 3, Table S5 and 
S6) relative to the observed record. Changes in the fre- 
quency and drought integral (sum of monthly drought 
SPEI values) for drought events (<0.02 percentile and > 
8 months) were also similar for climates with simulated 
changes in temperature or precipitation (Fig. 3C and D). 
Our results also show a sharp increase in the frequency of 
heat wave events (74%) and an increased likelihood of a 
heat wave occurring during an extreme drought event 
(mean increase of 113% Table S5) associated with simu- 
lated changes in temperature. 

Discussion 

Key attributes of climatic drought stress were assessed 
using a probabilistic approach that allowed us to describe 



the conditions capable of inducing drought damage and 
mortality in a consistent manner across a range of vegeta- 
tion types. Droughts coinciding with the documented 
die-off events reached minimum SPEI values of less than 
2% probability of occurrence within the observed record 
(Table 1). This finding suggests that regardless of regional 
climatic differences, tree populations among many species 
in Australian ecosystems tolerate at least 98% of the cli- 
matic conditions they experience and become vulnerable 
to drought stress events beyond this common climatic 
threshold. The occurrence of high maximum tempera- 
tures (MTI > 0.98 percentile in the vast majority of cases) 
and heat waves during all of the documented events rein- 
forces the importance of both water and heat stress acting 
in concert to bring about severe drought damage and 
die-off. Based on the climate profiles of the documented 
die-off events, we estimate that conditions that permit 
extreme droughts occur, on average, every 24 years (rang- 
ing from 14 to 62 years across sites). Robust and consis- 
tent estimates of drought occurrence provide an 
important benchmark for understanding the significance 
of drought in shaping vegetation patterns and modeling 
future changes in drought impacts. 

It is important to stress that these thresholds of drought 
intensity and high temperatures define an enhanced likeli- 
hood of tree die-off (above background rates) for a given 
site/region experiencing a similar drought history. Once 
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Figure 2. Projected changes in duration and 
intensity of future drought events relative to 
the historic climate for 15 tree die-off sites. (A) 
Drought duration versus intensity for observed 
(1891-2010) and future (2011-2050) datasets 
(using CSIRO Mk 3.5, A2 SRES scenario). 
Kernel density estimates of historic (solid blue 
line) and future (dashed red line) for (B) 
drought duration, (C) mean drought intensity, 
and (D) mean maximum temperature intensity 
(Mil). The mean maximum temperature 
intensity during each event is given by the size 
of the circle. The droughts associated with tree 
die-off events are shown as green squares and 
lowercase letters correspond to the seven sites 
shown in Fig. 1A-G. 
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drought events surpass these thresholds, the magnitude of 
the drought attributes (e.g., duration) and second-order 
factors such as soil type, topography, species, stand 
dynamics and disturbance history will affect the extent and 
pattern of die-off For example, it is likely that significant 
portions of the landscape (e.g., mid slope and valley floor) 
may only exhibit mild changes in vegetation structure and 
function, that is, reduced leaf area, while more sensitive 
sites (e.g., upper or exposed slopes) may sustain severe 
drought damage and tree die-off (Matusick et al. 2013). 
However, given the opportunistic and inherently biased 
nature of die-off observations (i.e., times when the presence 
or absence of die-off is often not reported) and lack of 
long-term monitoring in most ecosystems, the data used 
for this analysis cannot provide deterministic predictions of 
tree die-off for an extreme drought event without a more 
detailed appraisal of the sensitivity of the system through 
time and space. Nevertheless, this approach does provide 
insights into long-term and regional-level patterns in vege- 
tation responses to extremes in water deficit and tempera- 
ture and provides empirical estimate of likely return 
intervals for extreme droughts across a broad suite of eco- 
systems in Australia. 

The coincidence of high temperatures and heat waves 
with severe drought such as those documented here and 
elsewhere (Williams et al. 2013) is likely to escalate even in 
the absence of systematic declines in precipitation (Fig. 3). 
It is well established that summer rainfall and surface air 



temperature are negatively correlated at global and conti- 
nental scales (Lyon 2011; Mueller and Seneviratne 2012; 
Vautard et al. 2007). Thus, heat waves are common during 
drought, and the combined effect of reduced summer rain- 
fall and heat waves creates greater stress on plants than 
either factor alone (De Boeck et al. 2011). Average maxi- 
mum air temperatures were elevated during drought 
months (Fig. lA-G, Table 1), but more importantly, short 
periods of heat stress during summer were likely to have 
been pivotal in initiating tree die-off and mortality pro- 
cesses, when tree water and carbon status were already 
diminished (AUen et al. 2010; Galiano et al. 2011). During 
drought conditions, trees are more vulnerable to heat 
stress due to the lack of transpiration-mediated cooling of 
the foliage, enhanced heat loads on tissues, and increased 
light stress (Barua and Heckathorn 2006). Furthermore, 
increased evaporative demand will accelerate the decline in 
tree water status and promote dehydration of plant tissues. 

The lack of consistency in drought conditions preceding 
an observed event suggests that the antecedent meteorolog- 
ical conditions (assessed at 24 and 48 months) may 
provide limited prediction of the occurrence of die-off, 
suggesting additional biotic and abiotic factors not consid- 
ered here are likely to be important for vegetation sensitiv- 
ity to drought. For example, Fensham (1998) note that the 
persistence of wet periods over decadal and interdecadal 
periods can increase vegetation cover in northeastern 
Australian savannahs, thereby enhancing competition for 
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Table 2. Percentage change in projected or future (2011-2050) drought events relative to historic data (1891-2010) for all die-off sites in the 
analysis. Historic and future drought attributes are based on drought analysis of daily data from CSIRO Mk 3.5, A2 SRES scenario, medium sensi- 
tivity. Drought integral refers to the sum of monthly SPEI values during the event. Extreme events are defined as those with duration >8 months 
and where monthly SPEI and Mil values <0.02 and >0.98 percentiles, respectively. The return inten/al refers to the mean number of years 
between drought events. 
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Table 3. Details on 


the type of model used to fit the joint drought intensity and maximum temperature distribution for the observed and pro- 


jected or future (based on the CSIRO Mk 3.5, A2 SRES scenario) for each site. For the observed data, the probability of extreme conditions repre- 
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Figure 3. Future changes in temperature or precipitation have an 
approximately equal effect on drought events. Blue circles on all plots 
show the relationship between mean drought intensity and duration 
for drought events (SPEI <0.02 percentile and duration > 8 months 
for observed climate; 1961-2010). Red circles show simulated climate 
data perturbed by (A) temperature alone and (B) precipitation only 
(based on projections at 2050) for the 15 die-off sites (green squares 
denote drought associated with die-off events). The kernel density 
estimates of the drought integrals (sum of monthly SPEI values during 
drought) for all events based on obsen/ed (solid blue line) and 
perturbed climate data (red dashed line) for (C) temperature and (D) 
precipitation. Dashed reference lines (A-D) denote drought integral 
values at -10, -30, and -50. 



water and making these communities more vulnerable to 
die-off during subsequent drought events. The reverse may 
also be true, whereby extreme drought in the preceding 
years helps to reduce mortality for subsequent events and 
help to stabilize the system (Lloret et al. 2012). 

This study characterizes the joint probability density 
function of the two key determinants of climatic stress to 



provide a realistic assessment of current and future risks 
to tree survival, across their climatic niche and at the 
edges of their distribution (Fig. 1). By populating the bio- 
climatic space with observations of drought impacts and 
mortality, the limits of plant adaptation to drought can be 
quantified across a two-dimensional climatic space. While 
there are some similarities with this approach and other 
environmental envelope modeling involving associations 
among a suite of average climate parameters to evaluate 
species distributions, they are not directly comparable. An 
important element of our approach was to select climatic 
indices that have a direct functional role in limiting sur- 
vival rather than using a suite of climate variables that 
may only be correlated with species distribution. This 
means that our analysis is likely to be more robust in 
defining biologically relevant thresholds for episodic type 
die-off events. The characteristics of the copula distribu- 
tions themselves tell us how the likelihood of dry and hot 
conditions change as droughts become more extreme. The 
majority of these distributions conformed to the Gaussian 
or normal multivariate distribution with SPEI and MTI 
values showing differing levels of covariation. This tail 
dependence (Student-f copula), or lack thereof (Frank 
copula) in about one-third of the sites, implies that the 
co-occurrence of hot, dry conditions characteristic of 
extreme events occurs at proportionally higher likelihood, 
in the case of the Student-f distribution, or lower likeli- 
hood, in the case of the Frank distribution likelihoods, 
than would be predicted for their univariate distributions. 
For some ecosystems an increased or decreased coinci- 
dence of water deficit and high temperatures as part of a 
changing drought regime, will be a strong selective force 
for defining species adaptation to extreme drought events. 
We show that this multivariate approach is necessary for 
capturing the interactive effects of climatic stress and pro- 
viding a means to quantify changes in extreme drought 
conditions that have been derived from multisite analyses. 
(Fig. 1 and Table 3). This type of approach will help to 
resolve questions around the extent to which historic and 
recent changes in extreme conditions are symptomatic of 
a climatic shift and provide a means to assess changes in 
the probability of climatic thresholds likely to influence 
key issues around ecosystem dynamics, for example terres- 
trial carbon stores (Schwalm et al. 2010). 

The analysis establishes the predominant patterns in 
the observed drought regime within the duration versus 
intensity space and allowed us to assess changes in these 
attributes into the future (based on a moderate warm and 
relatively dry scenario) across different regions and vege- 
tation types (Fig. 2). In general, drought intensified under 
the projected future climate across all sites, yet on a site- 
by-site basis this intensification tended to be associated 
with increases in duration or intensity or both (Table 2). 
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Quantifying changes in the key attributes of drought 
events such as duration and intensity is critical for under- 
standing the mechanisms that underpin differences in 
responses that occur among individuals, species, and 
communities during severe droughts (Fensham and Fair- 
fax 2007; Koepke et al. 2010). Our current understanding 
of mortality mechanisms maintains that plant water and 
carbon balances are strongly linked to drought intensity 
by modulating soil and atmospheric water deficit across 
the soil-to-leaf hydraulic pathway and duration given the 
rate-limited nature of carbon assimilation and depletion 
(Mitchell et al. 2013b). Thus, the combination of inten- 
sity and duration of an "extreme" type drought, as 
described here, may have differential impacts on a partic- 
ular species, population, and/or cohort. While the 
approach presented here only describes the primary mete- 
orological drivers of droughts capable of inducing tree 
die-off, these data are crucial inputs for predicting site 
water balance in more detailed models that describe plant 
responses to drought events (McDowell et al. 2013). For 
example, shorter and more intense droughts are likely to 
kill trees via hydraulic failure, a process that has been 
widely attributed to episodic die-off events (Anderegg 
et al. 2012; Martinez- Vilaka and Pifiol 2002; Nardini 
et al. 2013). Alternatively, an increase in long-duration 
events at similar or reduced intensity may restrict carbon 
assimilation in species with relatively conservative stoma- 
tal regulation and potentially leading to carbohydrate 
starvation (McDowell 2011; Poyatos et al. 2013). Perhaps 
more importantly, prolonged droughts and elevated tem- 
peratures also increase the risk of biotic attack, while 
plant's defenses are compromised and/or conditions pro- 
mote greater activity and expansion of pest populations 
(Mattson and Haack 1987; Raffa et al. 2008). Changes in 
the frequency of extreme drought under the scenario pre- 
sented here (Table 2) and elsewhere (Mpelasoka et al. 
2008) may also reduce vegetation resilience through time 
if a complete recovery of plant vasculature, carbohydrate 
status and defensive mechanisms is not realized in the 
intervening years between drought events. A small num- 
ber of predicted droughts fell outside the margins of the 
observed record and are perhaps indicative of "mega- 
drought" conditions, characterized by higher intensities 
and longer durations than have ever been observed in the 
historic record (Fig. 2A). The climatic drivers of these 
mega-drought events under the projections used in this 
study originated from both increased temperature and 
altered circulation patterns that create prolonged declines 
in precipitation. If realized, these climate events may gen- 
erate unprecedented, extensive die-off that could induce 
long-term shifts in vegetation structure and function. 

Results from the sensitivity analysis of drought to 
changes in precipitation and temperature demonstrate that 



moderate increases in temperature alone have the potential 
to significantly lengthen drought events and increase the 
risk of extreme events even in the absence of reductions in 
precipitation. From a physiological perspective, elevated 
temperatures not only accelerate dehydration by increasing 
water loss but can impinge on plant carbon balance as 
well. Water-stressed plants with significantly reduced car- 
bon assimilation may deplete stored carbohydrates more 
rapidly when temperatures are elevated (Adams et al. 
2013). In turn, this may reduce plant fitness and impair 
recovery following drought, as carbohydrates are required 
for new growth during recovery from drought (Brodribb 
et al. 2010; Galiano et al. 2011). Changes in carbon metab- 
olism in response to elevated temperatures would have sig- 
nificant effects on ecosystem carbon budgets, by increasing 
respiration rates and the duration over which the ecosys- 
tem remains a potential source of atmospheric CO2 (Zhao 
and Running 2010; Keith et al. 2012). 

This study shows that the climatic drivers of episodic 
tree die-off across vegetation types and biomes in Australia 
can be characterized using a common drought surface and 
that tree die-off occurs beyond common thresholds in 
these indices. This study suggests that plant functioning 
adapts to variability in climate and the prevailing drought 
regime over a similar range of conditions. If this pattern 
were consistent globally, it would lend support to the 
concept that woody species from both wet and dry envi- 
ronments are similarly vulnerable to changes in water 
availability and temperature (Choat et al. 2012). The data 
presented in this study suggest that "top-down" 
approaches can yield important data on climatic thresholds 
of dominant species in their respective ecosystems, can be 
readily coupled to models of ecosystem dynamics from the 
local to global scale, and used to assess future mortality risk 
in forest systems. The attributes of drought events inducing 
less severe impacts such as changes in net primary produc- 
tivity and recruitment can also be assessed under this 
approach, thereby providing a means to characterize 
drought consequences on a range of ecosystem processes. 
Although the impact of extreme events such as drought 
may not always invoke a response in a plant community of 
similar extremity (Smith 2011), the nature of changes in 
drought conditions shown here may be rapid and severe 
enough to have significant and lasting impacts on species 
distributions and ecosystem dynamics. The phenotypic and 
genotypic plasticity of tree species, and associated stand- 
level responses to shifts in water availability, are the key 
determinants of future changes in community composition 
and distribution. If the thresholds associated with tree sur- 
vival are relatively fixed, as suggested by physiological evi- 
dence surrounding hydraulic failure (Martinez- Vilalta and 
Piiiol 2002), then the likelihood of drought events crossing 
these thresholds and inducing mortality wUl increase 
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significantly under future climate scenarios for many forest 
and woodland ecosystems globally. 
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Additional Supporting Information may be found in the 
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Figure SI. Drought history and details of the 2010-2011 
drought event near Jarrahdale in Australia's south west. 
Figure S2. Comparison of four GCMs and three scenarios 
for different sites across Australia. 

Table SI. Details of the 15 documented die-off sites used 
in this study (sorted by climate wetness; precipitation 
divided by potential evaporation). 

Table S2. Statistical parameter estimates from the joint 
distribution selection routine including; the log likelihood 
ratio, tau and degrees of freedom. 

Table S3. Percentage change in drought attributes 
between observed data perturbed by temperature change 
(based on CSIRO Mk 3.5, A2 SRES scenario, ~1.44°C 
change from 2010 climate) relative to observed data 
(1961-2010) for all documented die-off sites. 
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